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Abstract: We study the interaction between the AIMD (Additive Increase MultipHcative 
Decrease) congestion control and a bottleneck router with Drop Tail buffer. We consider 
the problem in the framework of deterministic hybrid models. First, we show that the 
hybrid model of the interaction between the AIMD congestion control and bottleneck router 
always converges to a cyclic behavior. We characterize the cycles. Necessary and sufficient 
conditions for the absence of multiple jumps of congestion window in the same cycle are 
obtained. Then, we propose an analytical framework for the optimal choice of the router 
buffer size. We formulate the problem of the optimal router buffer size as a multi-criteria 
optimization problem, in which the Lagrange function corresponds to a linear combination of 
the average goodput and the average delay in the queue. The solution to the optimization 
problem provides further evidence that the buffer size should be reduced in the presence 
of traffic aggregation. Our analytical results are confirmed by simulations performed with 
SimuHnk and the NS simulator. 

Key- words: TCP/IP modeling, router buffer sizing, deterministic hybrid model approach, 
multi-criteria optimization 
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Convergence et choix optimal de la taille du buffer du 
routeur pour le controle de congestion AIMD 

Resume : Nous etudions I'mteraction entre le controle de congestion AIMD (Additive 
Increase Multiplicative Decrease) et le routeur avec le buffer de type Drop Tail. Nous 
considerons ce probleme dans le cadre des modeles hybrides deterministes. D'abord, nous 
prouvons que le modele hybride de I'interaction entre la controle de congestion AIMD et 
le routeur de goulot d'etranglement converge toujours a un comportement cyclique. Nous 
caracterisons les cycles. Des conditions necessaires et sufRsantes pour I'absencc dcs sauts 
multiples de la fenStre de congestion dans le meme cycle sont obtenues. Puis, nous proposons 
un cadre analytique pour le choix optimal de la taille du buffer du routeur. Nous formulons le 
probleme du choix optimal de la taille du buffer du routeur comme probleme d'optimisation 
multi-critere, dans lequel la fonction de Lagrange correspond a une combinaison lineaire 
du taux moyen de transmission et le delai moyen dans le buffer. La solution au probleme 
d'optimisation fournit davantagc d'cvidcncc que la taille du buffer du routeur doit etre 
reduite en presence de I'agregation du trafic. Nos resultats analytiques sont confirmes par 
des simulations effectuees avec SimuHnk et le simulateur NS. 

Mots-cles : le modele de TCP/IP, choix optimal de la taille du buffer, modeles hybrides 
deterministes, I'optimisation multi-critere 
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1 Introduction 

Most traffic in the Internet is governed by TCP/IP (Transmission Control Protocol and 
Internet Protocol) [HHl]. Data packets of an Internet connection travel from a source node 
to a destination node via a series of routers. Some routers, particularly edge routers, expe- 
rience periods of congestion when packets spend a non-negligible time waiting in the router 
buffers to be transmitted over the next hop. TCP protocol tries to adjust the sending rate of 
a source to match the available bandwidth along the path. During the principle Congestion 
Avoidance phase the current TCP New Reno version uses AIMD (Additive Increase Multi- 
plicative Decrease) binary feedback congestion control scheme. In the absence of congestion 
signals from the network TCP increases congestion window linearly in time, and upon the 
reception of a congestion signal TCP reduces the congestion window by a multiplicative 
factor. Congestion signals can be either packet losses or ECN (Explicit Congestion Notifica- 
tions) [2^. At the present state of the Internet, nearly all congestion signals are generated 
by packet losses. Packets can be dropped either when the router buffer is full or when AQM 
(Active Queue Management) scheme is employed [lO]. Given an ambiguity in the choice of 
the AQM parameters [7l[T6], so far AQM is rarely used in practice. On the other hand, in 
the basic Drop Tail routers, the buffer size is the only one parameter to tune apart of the 
router capacity. In fact, the buffer size is one of few parameters of the TCP/IP network 
that can be managed by network operators. This makes the choice of the router buffer size 
a very important problem in the TCP/IP network design. 

The paper is composed of two principle parts. In the first part (Sections [2][5]) we analyze 
the interaction between the AIMD congestion control and the bottleneck router with Drop 
Tail buffer. This interaction can be adequately described by hybrid modeling approach. 
There are several hybrid models of the interaction between TCP and the bottleneck router 
[HElITSj. Here we analyze the model of [l3j. To our opinion, this model takes into account 
all essential details of TCP and at the same time leads to a tractable analysis. We show that 
the system always converges to a limiting behavior. In particular, we demonstrate that two 
different limiting regimes can coexist and the convergence to one or to the other depends on 
the initial conditions. Then, we provide necessary and sufficient conditions for the absence 
of subsequent packet losses. The absence of subsequent packet losses benefits the TCP 
performance as well as the quality of service for end users. We note that in [l^ there is no 
characterization of limiting regimes. Furthermore, in yjj only a sufficient condition for the 
absence of multiple jumps was obtained and the sufficient condition of [l^ is loose for some 
values of the decrease factor. 

In the second part of the paper (Sections [6]17|) we study the optimal choice of the buffer 
size in the bottleneck routers. There are some empirical rules for the choice of the router 
buffer size. The first proposed rule of thumb for the choice of the router buffer size was to 
choose the buffer size equal to the BDP (Bandwidth-Delay Product) of the outgoing link ^23j . 
This recommendation is based on very approximative considerations and it can be justified 
only when a router is saturated with a single long-lived TCP connection. The next apparent 
question to ask was how one should set the buffer size in the case of several competing TCP 
connections. In [5j it was observed that the utilization of a link improves very fast with 
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the increase of the buffer size until a certain threshold value. After that threshold value 
the further increase of the buffer size does not improve the link utilization but increases the 
queueing delay. Then, two contradictory guidelines for the choice of the buffer size have 
been proposed. In [17] a connection-proportional buffer size allocation is proposed, whereas 
in [3] it was suggested that the buffer size should be set to the BDP of the outgoing link 
divided by the square root of the number of TCP connections. A rationale for the former 
recommendation is that in order to avoid a high loss rate the buffer must accommodate 
at least few packets from each connection. And a rationale for the latter recommendation 
is based on the reduction of the synchronization of TCP connections when the number of 
connections increases. Then, [3l|T7] were followed by two works OITl] which try to reconcile 
these two contradictory approaches. In particular, the authors of [8] recommend to follow 
the rule of [3] for a relatively small number of long-lived connections and, when the number 
of long-lived bottlenecked connections is large, to switch to the connection-proportional 
allocation. One of the main conclusions of [TT] is that there are no clear criteria for the 
optimization of the buffer size. Then, the author of [11] proposed a general avenue for 
research on the router buffer sizing: "Find the link buffer size that accommodates both 
TCP and UDP traffic." We note that UDP (User Datagram Protocol) does not use any 
congestion control and reliable retransmission and it is mostly employed for delay sensitive 
applications such as Internet Telephony. We refer the interested reader to [M] and references 
therein for more information on the problem of optimal choice of buffer size. 

All the above mentioned works on the router buffer sizing are based on quite rough 
approximations and strictly speaking do not take into account the feedback nature of TCP 
protocol. Here we propose a mathematically solid framework to analyze the interaction 
of TCP with the finite buffer of an IP router. In particular, we state a criterion for the 
choice of the optimal buffer size in a mathematical form. Our optimization criterion can 
be considered as a mathematical formalization of the lingual criterion proposed in [llj . 
Furthermore, the Pareto set obtained for our model allows us to dimension the IP router 
buffer size to accommodate both data traffic and real time traffic. 

All proofs are provided in the Appendix. 



2 Mathematical model 

The window based binary feedback congestion control can be described by two functions 
f{w) and G(w). Function f{w) defines the increase profile of the congestion window and 
function G{w) represents the reduction of the congestion window upon the reception of 
congestion notification. Namely, in the absence of congestion notification the evolution of 
congestion window is described by the differential equation 

dw ^ f{w) 
dt T + x{t)/fi' 

where T is the two way propagation delay, x{t) is the amount of data in the bottleneck 
queue and fi is the capacity of the bottleneck router. We note that T + x{t)/fi corresponds 
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to the Round Trip Time (RTT) when the amount of the enqueued data in the bottleneck 
router buffer is x{t) at time moment t. Thus, function f{w) determines the increase of the 
congestion window per one Round Trip Time. The sending rate X{t) of the window based 
congestion control is given by 

We would like to emphasize that here the time parameter t corresponds to the local time 
observed at the router. 

We study a Drop Tail buffer with size B. \{x(t) < B, the congestion window w increases 
according to ([T]). When x reaches B at time t* , i.e. x{t*) = B, the buffer starts to overflow. 
The overflow of the buffer will be noticed by the sender only after the time delay S = T+B/fi. 
Upon the reception of the congestion signal at time t* + S, the congestion window is reduced 
according to 

w{t* +S + 0) = G{w{t* +5-0)). (3) 

As we shall see below, w (resp.. A) can represent either a congestion window (resp., sending 
rate) for a single TCP connection or a total window (resp., total rate) of several TCP 
connections. 

Consider n long-lived AIMD TCP connections that share a bottleneck router. Denote by 
Wi(i) the instantaneous congestion window of connection « = 1, ...,n at time t e [0, oo). In 
the case of the AIMD congestion control, if x < B the evolution of the congestion window Wi 
is given by differential equation |[T|) with f{w) = = const. If we restrict ourselves to the 
symmetric case Ti = T — const and rrii = mo = const, the sum of all congestion windows 
w(t) = X]r=i^^(^) ^^^^ satisfies differential equation |(T|) with f{w) — m, where m = nmo. 
Namely, we have 

dw m 

'dt ^ T + x(t)/fi' ^ ' 

dx _ ( X{t) - fi, if < x{t) < B, or x{t) = and \{t) > or x{t) = B and X{t) < ^; 
dt \ otherwise, 

(5) 

where X{t) is given by And iix{t*) — B at some time moment t* , the congestion window 
is decreased multipHcatively after the information propagation delay 6 = T+B/^as follows: 

w{t* +6 + 0) ^ P'''w{t* +5 -0), (6) 

and consequently, G{w) = (3^w for the AIMD case. Usually, k = 1, but sometimes it is 
necessary to send several congestion signals in order to reduce the sending rate below the 
transmission capacity of the bottleneck router. 

Since we consider the case of equal propagation delays, the synchronization phenomenon 
takes place [TOj, and consequently, the total sending rate is also reduced by the factor 
For instance, in TCP New Reno version the reduction factor (3 is equal to one half. 
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Let us make the change of time scale according to 

A dt 



ds 



and the change of variables: 



A , A , 

V = w/m, y = x/m. 



The new time s can be viewed as a counter for Round Trip Times. Now the dynamics of 
the system between the jumps is described by equations 

^ _ / ^(0 ~ yit) ~ 9i if < y{t) < b, or y{t) = and v{t) > q, or y{t) = b and v{t) < q + b; 
ds \ otherwise, 

(8) 

where q = fiT/m is the maximal number of packets that can be fit in the pipe, in other 
words Bandwidth-Delay Product (BDP) in packets, and b = B/m is the maximal number 
of packets that can be fit in the router buffer. Let s* be the moment in the new time scale 
when component y reaches value b. Then, equation ^ is transformed to 

i;(s* + l + 0) = /?S(s* + l-0), (9) 

where k = m\n{i : [3'^v{s* + 1-0) <b + q}. 

Remark 1 Because of the delay in the information propagation, the congestion window is 
reduced after the delay 5 — T + B/fi in the original time scale, or, equivalently, after 1 time 
unit in the new time scale s. The value of k is such that, after sending k congestion signals, 
the amount of data x {and y ) starts to decrease. 



3 Convergence of the system trajectories 

The dynamics is defined by three parameters /3, g, and 6, and the system trajectory remains 
in the region ^ = {0 <y <b, w > 0}, provided the initial condition is there. 

Suppose a trajectory starts at s = from initial condition j/o = b, (3{q + b) < vq < b + q^ 
and s* is the first moment when y{s*) — b. Let vi — v{s* + 1 + 0). We introduce mapping 

(p such that vi = (p{vo)- Consider the iterations u^+i = ^{vi), i = 0, 1, .... 
Theorem 1 There exists limi^ooVi ~ ^(wo) with 

^ ' \ V"2, i/v e {d,q + b), ^ ^ 

for some constant d. In particular, one of the above intervals can be empty. 



^Initial conditions outside the region [/3(q + b),q + b) are of no interest because, after the very first 
(multiple) jump we have v{s* + 1 + 0) £ [/3(g + b),q + b). 
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Definition 1 Suppose the trajectory starting at s — from initial condition yo = b, vq < 
b + q reaches the same point, for the first time, at some time moment S > 1. Then this 
finite trajectory is called a cycle. A cycle with component y remaining zero for a positive 
time interval is called clipped (see FigureU^. If a cycle touches the axis y — only at a 
single point, we call such cycle critical (see Figure\^. 

Corollary 1 (from Theorem [l| Any cycle has a single time moment, when a (multiple) 
jump occurs. 

The number k of instant jumps of component v is called a cycle order. We call such 
cycles /c-cycles for brevity. If one of the intervals in ifTOl) is empty then only a single cycle 
exists (Figure [1]). Otherwise, two cycles exist simultaneously (Figure [S]); their orders are 
two subsequent positive integers. According to Theorem [1] which cycle is realized depends 
on the initial conditions. 



4 Properties of cycles 

In this section, we characterize the shape of cycles. In other words, for given parameters 
P, q and b, we would like to know if the limit cycles of the system trajectories are clipped 
or unchpped and what orders the cycles have. For fixed values of P and q, we define the 
following quantities: 

A^ = min|z>l: < g| ; (11) 



D^Hl-(3^) + ^^; (12) 



C = -ln(l-/3^)-/3^; (13) 
9k is the single positive solution to equation 

9 3^9 3^ 

1 - 6^"*= 1 - e^^fc 

Then, we define the set of quantities which do not depend on q: 
Tk is the single positive solution to equation 

l-e~^, fc = 2,3, ... (16) 



1-/3 



+ 
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It is convenient to put ti,A* and ql equal to +00. Finally, in case q < D one has to solve 
equation 

e-"^ + r-l^ I3^{q + r + l)-q. (19) 
It has no more than two positive solutions r<f which define 

b = e-^ + r-l; b^e-'^ + f-l, (20) 

Note that b<b. If g < g^r+i then q < D and b > A*j^_^^ — q. 

We note that all the above defined quantities do not depend on b. Thus, from now on 
we assume that /? and q are fixed and we are going to describe what kind of cycles exist for 
different values of b. In other words, we study what effect the router buffer size has on the 
limiting behavior of TCP/IP. There are three cases: 



Case A*j^,^ < q 



If 6 e 

13" 



N+1 



-P n 

iZTpo — -roo for generality. 



then only the cycle of order N exists. In case TV = 1, we put 



- q,A 

exist simultaneously. For b £ (^A*pf — <?, - 



N-1 



two cycles, of orders N and iV — 1 



, there exists only a single cycle of order 



— 1. And so on; for b > A2 — q, only 1-cycle exists (see Figure fT3l) . 
The A^-cycle is cHpped for b G [0, 60, w)- Cycles of lower orders are uncHpped for all values 
of b, if they exist. 

The A^-cycle touches the f-axis at a single point iff 6 = 60, Af- Thus, if 6 = 60, Af there 
exists a critical A^-cycle. No critical cycles of lower orders exist. 

Example 1 Let us illustrate this with a numerical example. If we take q = 0.9 and (3 = 1/2 
then N = 2, Al = 1.4965, Al = 0.3910. If b € [0,0.1] we have only 2-cycles; if b € 
(0.1,0.5965] we have 1-cycles and 2-cycles {see Figure\^; and if b > 0.5965 we have only 
1-cycles. For each b < 60,2 = 0.0617, there exists only a clipped 2-cycle [see Figure\^. As 
one can see on Figure\^ when b — 60,2 — 0.0617, the 2-cycle becomes critical. All figures 
for this example have been plotted with MATLAB Simulink. 



Case q < g^+i- 

If 6 G [0,^), then only the TV-cycle exists. If 6 G [b,A*pf^^ — q], then two cycles of orders 

A'' and A'' + 1 exist simultaneously. For b e (^jv+i — q, i_pN-i — q], again, only the A''-cycle 
exists. 

The A^-cycle is clipped for b e [0, 6o,Jv); the (N + l)-cycle is clipped for b G [&, 6o,Af+i)- 
These cycles become critical at b = 60, Af and b = 6o,Af+i; respectively. Cycles of lower orders 
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Figure 1: Clipped 2-cycle. Case A2 < q. 



Figure 2: Critical 2-cycle. Case A2 < q. 
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Figure 3: Co-existence of 1-cycle and 2-cycle. Case < q. 



are undipped for all values of b, if they exist. If > 1 then, similarly to the case ^^^i < q, 
the order of the cycle decreases as b increases above ^_^jv-i — q (see Figure [13]) . 
Case g^^i <q< ^jv+i- 

If C < A*j^_^j^ — 5 □ and q < D, then everything is similar to the case q < q'^^i- The 
difference is that the (iV -t- l)-cycle is cHpped and cannot be critical; it exists simultaneously 



only the A^-cycle exists. The latter 



with the TV-cycle for b € [b,b]. If b e (b, ^3"-^ 
interval is non-empty. 

If C > A*j^_^^^ — q or D < q, then everything is exactly as in case A*j^^^ < q. 



Actually C cannot be equal to A*^^-^ 
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5 Conditions for the absence of multiple jumps 

The regime with multiple jumps is not desirable. The multiple jump corresponds to the 
lost of more than one packet in a single congestion window. Subsequent packet losses can 
force TCP to switch from the Congestion Avoidance TCP phase to the Slow Start phase and 
lead to lengthy timeouts. Furthermore, the absence of subsequent packet losses is beneficial 
not only for the TCP performance but also for the quality of service provided to the end 
users. In the next theorem we provide necessary and sufficient conditions for the absence of 
multiple jumps, namely, we characterize all possible cases when only a single cycle of order 
1 exists. 

Theorem 2 The following mutually exclusive conditions fully characterise all possible cases 
when only a single cycle of order 1 exists: 

(a) > Q ind b + q > Ai,; 

(b) A2 < q (b can be arbitrary); 

(c) ^<q<q*2 and b ^ [b, A^ - q]; 

(d) max I j4a, q^^ < q < A^ - C , q < D and b <^ [6,6]; 

(e) max | jr/ji < 9 — ^2 ^ 9 > can be arbitrary); 

(f) max|Y3g, gji ^2 ^ c| < 9 5; ^2 can be arbitrary). 

In the following corollary we provide a simple sufficient condition for absence of multiple 
jumps. 

Corollary 2 Condition b + q > A2 is sufficient for the absence of cycles of orders k > 1. 
( See Figure [HI and Corollary O ) 

Recall that A^ depends only on (3. In particular, if /3 = 1/2, A^ = 1.4965 . 

We would like to note that the above sufficient condition is tighter than the sufficient 
condition for the absence of multiple jumps provided in [l3j: b -\- q > 2/3/(1 — fi). To 
compare these two conditions, we plot A2{I3) and 2/3/(1 — /3) in Figure [Hand the difference 
2/3/(1 — /3) — A^ifi) in FigurelH Strictly speaking we have 

Proposition 1 The difference 5 ^ — ^2 always positive and lim^^i 5 ~ +00. 

Nevertheless, the simple sufficient condition of [13] appears to be quite good except for 
values of (3 that are too close to one. 
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6 Pareto set for optimal buffer sizing 

Let us study what effect has the choice of the buffer size on the performance of TCP. In 
particular, we are interested in optimal buffer sizing. Towards this goal, let us formulate the 
performance criteria. On one hand, we are interested to obtain as large goodput as possible. 
That is, we are interested to maximize the average goodput 

1 /•* 

g = hm - / g{s)ds, 

where the instantaneous goodput g{t) is defined by 

_ r A(i), if x{t) < B, 
' ~ \ if x{t) = B. 

On the other hand, we are interested to make the delay of data in the buffer as small as 
possible. That is, we are also interested to minimize the average amount of data in the 
buffer 

t 



1 f 

lim — / x{s)ds. 

:^00 t Jq 



t— >oo t 

Clearly, these two goals are contradictory. In fact, here we face a typical example of multi- 
criteria optimization. A standard approach to it is to consider the optimization of one 
criterion under constraints for the other criteria (see e.g., |JJJ). Namely, we would like to 
maximize the goodput given that the average amount of data in the buffer does not exceed 
a certain value 

max{5 : x < x*}. (21) 
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Or we would like to minimize the average delay given that the average goodput is not less 
than a certain value 

fmm{x:g>g^}. (22) 

The solution to the above constrained optimization problems can be obtained from the 
Pareto set. As is known, see e.g. [19], the Pareto set can be constructed by solving the 
optimization problem 

1 



lim 

t^oo t 



cig{s) — C2x{s)ds 



(23) 



To be more precise, the Pareto Set is formed by the pairs of objectives {g, x) that solve ((23|) 
for different (ci, C2) G R\. An example of Pareto set is given in Figure [6l Each point of the 
Pareto set corresponds to a solution of optimization problem l(23|) for some choice of ci and 
C2. Once we obtain the Pareto set, it is very easy to deduce solution of problems l(2T|) and 
l(22|) . For instance, if one wants that the utilization of the bottleneck router will be not less 
than, say, 95%, one has to be ready to accept the delays that are equal or greater than x,. 



1^ 

0.95|a. 







//////////// 



X , 



Figure 6: Pareto set. 








D 



q ^ 

Figure 7: Phases of the clipped cycle. 



All three optimization problems (|2T|) . l|22|) and l[23|) can be regarded as mathematical 
formulation of the lingual criterion "find the link buffer size that accommodates both TCP 
and UDP traffic" given in [ll|. Since UDP traffic does not contribute much in terms of the 
load, for the design of IP routers one can use for instance optimization problem l(2T|) where 
the delay constraint is imposed by the UDP traffic. 

We note that here we deal with the optimal impulse control problem of a deterministic 
system with long-run average optimality criterion. To the best of our knowledge there are 
no available results on such type of problems in the literature. In principle, the control 
policy in our model can depend on the current values of x and A. In practice, however, all 
currently implemented buffer management schemes (e.g., AQM, Drop Tail) send congestion 
signals based only on the state of the buffer. Thus, we also limit ourselves to the case when 
the control depends only on the amount of data in the buffer. Furthermore, we restrict the 
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control action only to the choice of the buffer size. Thus, the control signal is only sent at 
the moment when the buffer gets full. 

The following theorem provides expressions for the average sending rate, goodput and 
queue size under condition q > A2, which guarantees the absence of multiple jumps for any 
value of the buffer size. Remember that A2 depends only on /? (see (fT6 |) . lfT7|) ). In particular, 
the expressions allow us to plot the Pareto set parameterized by the buffer size. 

Theorem 3 Let the condition ^T/m > A2 be satisfied. Then, for B e [0, m6o,i] the average 
sending rate, goodput and buffer occupancy are given by 



m{l- 0^) 



2T, 



cycle 



T, 



- I H SCD 

m 



mT 



cycle 
1 



1 + ^ + ScD 

m 



TO 



TO 



cycle 



Sae 



TO 



Sae 



VAB{s)ds 



yAB{s)ds 

Sen 



SCE 



ycr,{u)du 



ycD{u)du 
B{plT + B) 



respectively, where Tcyde is the cycle duration given by 



B to 



Tcycle = (1 - + — + Scd)T +- + 



Sae 



yAB{s)ds 



ScE 



ycD{u)du 



with 



ycD{u) = e " + (u- 1), 



VAsis) + + pScD]e-' + - 1) + (i{ScD + 1) - (1 - /?) — , 

TO 771 TO 

where Sen and Sab are the solutions of the equations 



-Sc 



" + S, 



CD 



B 

TO ' 



- - fiScD + {I - + —) 

TO TO 



For B e (to6o,i,oo), we have 



i-^^^ + Sab + pScD - (1 - /3)(1 + ^) = 0. 

TO 



2T, 



TO 1+ P 



cycle 
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9 = fJ; 

2 / psi 



where 



with 



mT / y{s)ds -\ ( / y'^{s)d 

Jo \Jo 

Tcycie = T{si + 1) + — ( / y{s)ds + — 



B{tiT + B) 



vis) = 



1 H vq 



e-' + (s -l)+vo- 



where vq and si are defined by l\28^ and l\29\i with k — 1. 



Example 2 Let us illustrate the Pareto set for a benchmark example of the TCP/IP network 
created with the help of NS-2 simulator [I^ . The network consists of a single bottleneck link 
of capacity /i = \QMhps which is shared by n long-lived TCP connections. The propagation 
delay for each connection is T = 0.24s and (3 — 1/2. The packet size is AOOObits. Thus, 
we have that tuq = AOOObits as well. In Figurel^ we plot the Pareto set for n = 10 (and 
m = nniQ — 40, 000^ using the formulae of Theorem\M and measurements obtained from NS 
simulations. As one can see, two curves match well. In Figure{^ again using the formulae 
of Theorem\^ we plot the average goodput and the average sending rate as functions of the 
buffer size for n = 60. 

We note that g < always, but the average sending rate A can exceed the router capacity 
/X (see Figure [9]). Nevertheless, as the next Proposition [2] states, the difference between the 
average sending rate and the router capacity goes to zero as B increases. In particular, this 
means that when the Drop Tail router is used, the rate of lost (and then retransmitted) 
information eventually diminishes to zero as the buffer size increases. 



Proposition 2 When B ^ oo, the difference A = A — /i approaches zero from above. 



7 Minimal buffer size for the full system utilization 

In the case of multiple TCP connections competing for resource of the bottleneck router we 
have m = nrriQ. Here n is the number of competing TCP connections. Let us study how the 
minimal buffer size for the full system utilization, BQ pf, depends on n or, equivalently, on 
TO. Bo, AT is the buffer size corresponding to scenario when the Pareto set touches the level 
/z (see Figure [6]). It corresponds also to the critical cycle of minimal order. 

Proposition 3 (a) For a fixed N , the value of i?o,JV = "mbo^N decreases as to increases. 
(6) The value of -Bq.tv increases as N increases. 
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NS-2 simulations - 

Theoretic curve 



average occupancy (packets) 



^— average sending rate 
■ - average goodput 



1000 1500 2000 2500 

buffer lengtfi, B (packets) 



Figure 8: Pareto set: Numeric calculations Figure 9: Non-monotonicity of the average 
and NS-2 simulations. sending rate. 

Corollary 3 The buffer size -Bq.tv of the minimal order critical cycle is a piece-wise differ- 
entiable function of m, decreasing on the intervals [m.i,mi+i); 

lim Bo^Nim) < -Bo.iv(m^+i), ^ = 0, 1,2, . . . 

m — ^rrii^i —0 

Here rrii = /iT(l — /3')//3*; the value of N equals i + 1 on the interval [mi, mi+i) [see ). 

Moreover, lim„j^™„_o Sq^w = 0, lim„i^„„_o = 0, lim„^o+ -So,i = t^T[l - P)/l3, 

limAT^oo nT-N-iBa^Ni'mN-i) — 0.5(/iT(l — /3))^ and hence lim„j^oo Bq n = 0. 

Example [2](cntd.) In Figure fiOl we plot the buffer size Bq^n of the minimal order critical 
cycle and the curve f{m) = (1 - l3f{nTf/{2m) for fiT = 2.4 x lO'^bits (GOOpackets) . The 
curve f(m) indeed approaches fast the local maxima of Bo pf as m increases. In Figure [771 
we make a zoom on the interval with smaller values of m. As one can see, when m goes to 
zero, the value of Bq n approaches SOOpackets, which is the BDP in this network example. 

We note that by Corollary [3] for small values of m the minimal buffer size for the full 
system utilization is approximately equal to jjT , BDP of the bottleneck link. This is in 
agreement with the empirical conclusion of [23]. In [3] the authors suggested that the mini- 
mal buffer size for the full system utilization should decrease as (/iT) / ^/n as the number of 
connections n increases. We note that the authors of [3| have assumed that the competing 
TCP connections are not synchronized. That is, only a single connection reduces its con- 
gestion window when the buffer becomes full. In our model we assume full synchronization 
of competing TCP connections. Namely, when the buffer is full, all connections simultane- 
ously reduce their congestion windows. We expect that the situation in real networks is in 
between these two extremes. And thus, the model of [3] provides an upper bound and our 
model provides a lower bound. Furthermore, it was beheved previously that if the compet- 
ing TCP connections are synchronized, one has to provide BDP of buffering to guarantee 
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the full system utilization. From Figure [TT] one can see that the minimal buffer requirement 
decreases with increasing m (or, equivalently, with increasing n) even in the case of complete 
synchronization. Finally, we would Hke to mention that the value of Bq.n is non-monotonous 
with respect to m, even though it eventually decreases to zero (see Figure [T0|) . Curiously 
enough, the experiments of [24] with the router, running FreeBSD dummynet software, have 
also shown the non-monotonous behavior of the minimal buffer requirement in the case of 
synchronized connections (see Figure 1 in [23]). 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 10 20 30 40 50 60 

" xio* n 



Figure 10: The minimal buffer size for the full Figure 11: The minimal buffer size for the full 
system utilization. system utilization (zoom). 



8 Conclusions 

In this paper we have studied the interaction between AIMD Congestion Control and a 
bottleneck router with Drop Tail buffer. We have used the hybrid modeling approach. It is 
demonstrated that the system always converges to a cyclic behavior. The limit cycles have 
been fully characterized. In particular, we have obtained necessary and sufficient condition 
for the absence of cycles with multiple jumps and a simple but tight sufficient condition. 
Then, we have formulated the problem of choosing the buffer size of routers in the Internet as 
a multi-criteria optimization problem. In agreement with previous works, our model suggests 
that as the number of long-lived TCP connections sharing the common link increases, the 
minimal buffer size required to achieve full link utilization decreases. However, in the case 
of synchronized connections, the decrease is not monotonous and slower than the inverse 
of the square root of the number of connections. The Pareto set obtained with the help of 
our model allows us to evaluate the IP router buffer size in order to accommodate real time 
traffic as well as data traffic. The simulations carried out with the help of Simulink and NS 
Simulator confirm the qualitative insights drawn from our model. Application of the same 
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framework to other congestion control mechanisms, such as MIMD, HighSpeed TCP, TCP 
Westwood appears to be a fruitful direction for future research. 



Appendix 

Undipped cycles. 

In this and the next subsection, we ignore the requirement that y > 0. Thus dynamics 
is described by equations 

dv 1 . 

ds ~ ^' 

( v-y-q, if y < &, or , . 

^ = 1 y = 6 and w < A; ^ ' 

[ otherwise, 

where 

A^b + q. 

The jumps occur according to (7) as before. 

Definition 2 Let yo — b and vq < A be the initial conditions. A piece of trajectory on the 
time interval [0, s* + 1 + 0] is called a pseudo-cycle of order k (see (7)). If v{s* + 1 + 0) = wq 
then the pseudo-cycle is called a k-cycle. 

Later, it will be shown that if a cHpped k-cyde exists then the undipped fc-cycle exists, 
too (Corollary [6l). Clearly, (|24l) has a single solution 



v{s) =vo + s; 

y{s) ^ {1 + q + yo- t;o)e"^ + s - 1 + vq - q. 



Theorem 4 An (undipped) k-cycle exists iff 



(25) 



(26) 



whe 



Al = I 1-/3" ' VK > I, ^27) 



oo, if k — 1 

and, for k > 1, Tk is the single positive solution to (7 



Proof. Obviously, parameters of a fc-cycle, vo and time interval si can be found from 
equations 

y(si) = 6; p''v{si + l) = vo, 
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which are equivalent to 



1 



13^ 
si 



1 -X- A ^"(''i+i) 

I ^ l^pk 



(28) 
(29) 



A fc-cycle exists iff ((29|) has a positive solution and vq given by ((28|) satisfies inequality 
> PA. (Otherwise, if vq < [3 A, there is no need to reduce v so many times.) Equation 
l(29|) has a positive solution iff 



d 



l + A 



1-/3'' 



< 1 



s=0 
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Lemma 1 If vq £ [(3A,A) then, starting from vq, yo — b, the next instant series of K + 1 
jumps results in the value v < A. Hence the order of any cycle cannot exceed K + 1 (and 
clearly cannot be smaller than K). 

Proof. Suppose vq — (3A. Then, after the next instant series oi K + \ jumps, the value 
V is not smaller than v. To put it differently, 

v<l3^+^[l3A + s+l\, (32) 

where s solves equation 

(1 + A - /3A)e"" +s-l + (3A = A 

If we substitute 



into l(33|) we obtain, using equality A = "'" 



l + e-^ = S:i4-l + e->^-l + e->0. 



1 + (1 " f3)A (3{2 + S) 1 + S 



When s increases from zero, the lefthand side of l|33|l initially decreases from zero and 
increases thereafter. Hence s > s and l(32|) implies 



V < p'^+^lpA + s + 1] < /3^+i[/3^ + -(3A- 1 + 1] = A. 



Lemma 2 Suppose j3 G (0, 1) is fixed and consider function 

fk{A)^{l-l3'^)A-(3^-\s, + l), (34) 

where si solves f^ff)) . The domain of f is given by hSOjl . Then 

(a) > 0; 

(b) fi{A)<Ofor allA>j^; 

(c) \fk > 1 equation fk{A) = has a single finite solution A'^ given by ^27]) : A'l decreases 
as k increases. 

(d) Vfc >1AI> Y^^; Vfc >2,Al< 
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Proof , (a) According to the rule of implicit differentiation, applied to equation 



we have 

dsi (1-6-^1)^(1-/9'') 



dA e-''isi(l - I3'') - (1 - 6-^1)^/3'' - (1 - - e-«i)' 

The denominator equals 

-(1 - e-"' - sie-"') - /?'=(sie-^i - e'^i + e'^^^ 

< (i^ie"^ - e-2«i _ sie-2"i) - (1 - e""' - sie^'O 
= (1 - e-"i - sie-"i)(/3e-''i - 1) < 0; 

hence ^ > for si > 0. 
Now 

dJk{A) ^ _ ak^ _ nk-l<tX 
dA P ) P dA 

_ (1 - I3^)[{1 - f3''){sie-'^ - 1 + e-"0 + p^-^jl - er'^ fjl - /?)] 

^ sie-''i(l - 13^) + (1 - e-'*i)(/3''e-"i - 1) ' 

The denominator is negative (see above). The nominator does not exceed 

(1 - /3'=)(1 - /3)[(sie"^i - 1 + e-'^)[l + fi'^-^) + /3'=-i(l - e-'^f] 

= (1 - /3'=)(1 - /3)[(sie-^i - 1 + + /j'^-ie-^^si - 1 + e"'')]- 

The both terms in the latter square bracket are negative for si > 0. Hence '''^''J^^^ > 0. 
(b) It is sufficient to prove that 

si > S = A{1- I3)-1, 

where si solves l[29|l at fc = 1. 

Case 5 < is trivial, thus assume that S > 0. Let us substitute S into the both sides of 
l(29|) and estimate the difference: 

S -, — ^ ^ -.— ^2 



1 I 4 /3(g+l) 1 I S+l ' 5' + 2 

^ 1-/3 1-/3 1-/3 

because function (5* + 2)e~^ decreases from 2 at 5* = 0. To complete this part of the proof, 
it is sufficient to notice that, on the interval 
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the righthand side of l|29p is smaller than the lefthandside iff S" < si. 

(c) The first part is obvious: A^. is given by l(27|) . provided equation (fTGl has a single 
positive solution. The latter statement folows from the fact that function 

ff(r) = (l-e--)(l + a(r + l))/r 

decreases to lim,-_>oo giT) = a, starting from limT-^o = 1 + a. Here 

A (3''-^ - 13'' 



1-/3'= 

Indeed, 



(35) 



dg _ €^''[1 + a + T{l + a + ar)] - {1 + a) ^ 
dr 

in case a < 1, and 

Now, look what happens as k increases. Obviously, functions f^rpk — i_pk ■ ^ and 

a = jzipkijj — 1) (see l|35p ) decrease. According to (f27|) it remains to prove that given by 
ifTG]) increases with a. We rewrite (fTGl as (1 + a{T + 1))(1 — e^'^) — r = 0. Hence 

dn ^ (Tfc + l)(l-e-"'') ^ (Tfc + l)^(l-e-"'°)^ 

da a{l - e--^") + (1 + a{Tk + l))e-^fc - 1 h{Tk) 

where Mt) = -2 + 3e"^ - e^^r _^ ^g-r ^2^-t_ ^laNe substituted a = ./"I^+^r"'; 
We intend to prove that 

^ = -2e-" + 2e-2- + re-^-r2e-^<0 (37) 

UT 

when r > 0. Clearly (l37| holds for r > 1. 
Suppose r e (0, 1). Then 

d'^h 

— = 3e-^ - 4e-2^ - 3re-^ + r^e"^ < 36"^ - 4e-2^ - 2re-^ = £"^(3 - 46"^ - 2r). 



Expression in the brackets has a negative maximum at t = In 2. Therefore, < and 
< 0. Finally, /i(r) < for all r > 0, because /i(0) = 0. 

(d) To estimate A*^. from below, we use statement (a): it is sufficient to establish that 
fk (iS^) < 0, ie Si + 1 > ^ ^, ,.-1^' ' 3/=(s+i T < 1 - e-^ for S = - 1. 

(The argument is similar to (b).) But 



1+, 3fc-l- i_^fc 



- 1 + e"^ = - 1 + e"'' = e"'' - < 



1 §1 1-/3" 5 + 1 1 + 5 
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because function e '^(1 + S*) decreases from 1 at S" = 0. 

Finally, in case k > 2, suppose A^. > i_pk^2 ■ Then for parameters values /? and A S 



A*f.j we have that ((30|) holds for k~2, k~ 1, and k and simultaneously A- (A) < 0, 
fk~i{A) < 0, fk-2{A) < 0: see (a) and (c). According to the beginning of the proof of 
Theorem 1, cycles of orders k, k — 1, and k — 2 exist which contradicts Lemma [TJ H 

Now we can easily finish the proof of Theorem IH Suppose a /c-cycle exists. Then, ac- 
cording to ([3Q1) . a > jzTgk- Lemma El guarantees that 



Al > 



1-/5 



> 



1-/3'= 



and, as was mentioned earlier, inequality vq > (3 A must be valid (see (|28l) ). which is equiv- 
alent to A < Al- Finally, if ((2 
'^o > P^', hence a fc-cycle exists. 



alent to A < Al. Finally, if ((26|) holds then ([29l) has a positive solution (see ((30l) ) and 



Case A 



N+i > <i- Case A*f^^^ < q: 



5-cycle exists | 



4- cycle 
exists 



3-pycle exists 



■7 



2-cycJe exists 




1-cycle exists 



A* - (J Al-q Al-q A*-q 

Figure 13: Existence of uncHpped cycles; = 4. 



Remember that A = b + q. Thus, if q is fixed and b increases from 0, undipped cycles 
have orders N (see ifTT]) ) and, possibly, TV + 1, if A'^^-^ — q > 0. Later, as b increases, the 
order of cycles decreases according to Fig. [131 

Stability of undipped cycles. 

We intend to study the mapping Lp introduced just before Theorem [TJ Since we study 
only undipped cycles, this map is a little different and will be denoted (f. But firstly we 
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concentrate on a different mapping: 



A 



defined for vq G [I3A, A] under a fixed fc > 1. Here s* = if = A; in case vq < A, s* > 
is the first moment when y{s*) = b starting from y{0) = b, v{0) — vq. 



Lemma 3 



dvo 



< (3^ and hence ^'^ is a contraction. Function is decreasing. 
Proof. Assuming that vq < A, s* is a single positive solution to equation 

(l+A-wo)(l-e-")-s = 0, 



(38) 



hence 



and 



Finally, 



ds* 



1 - e- 



1 - e- 



(1 + A - i;o)e-** -1 



avo dvo s*e 



1 + e- 



s* -1 + e- 



s*e ^ + e 



1 e" 



' - 1 
- 2s* 



< 0. 



1-e- 



> 0, 



because the nominator increases, starting from at s* = 0. Therefore > {—P'^). 



iff d < (3 A, where d is a solution to $ [d) — A. Here 



Lemma 4 (a) A [^k+i^ 

gK-l A I— II 

and below, ^_pK-i = <x if K = 1; K is defined by ^31]) . 

In this case, Vwq G [f^A,A), the mapping (p{vo) coincides with ^^{vq). 
(b) If A G ^ i-pi< ' ^*K+i following statements hold: 

(a) Vvo e [PA,d], <f{vo) = <f''+Hvo) e [PA,d]; 
((3) Wvo G id,A), 0ivo) = <i>^K) e id,A). 
See Fig\14\ (Note that, according to the definition of K , A E {^-^ 

to Lemma[M ^k+i ^ (t^^' 1-/3^-1 ■) 



; according 



Proof, (a) According to the definition, d — -0^ ^ s* — 1, where s* solves l(38|) under vq — d. 
\id = I3A then 

{l + A-(3A){\-e-'') = s*; 



— - s* - 1 



PA, 
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or equivalently 

^ — 1-0K+1 , 

To put it diflerently, we have A = A'^_^_-^ if d ~ l3A. 

It remains to prove that d — (3 A = — s* — 1 — (3 A is a decreasing function of A. Since 
s* satisfies equation 

(^l + A-^+s* + l^ (l-e-^*)-s*=0, 

ds^ _ (l-/?^)(l-e-^*)^ 
rf^ ~ /?^e-«*(e-«* + s* - 1) 

and 

d(d - PA) _ \_ ds* s*e-''' ^e-'" - \ + e"^")^ - d^+^er"' [e.-'" + s* - 1) 

'dA " W^'dA^ /3^e-«'(e^''* + s* - 1) 

The denominator is obviously positive for s* > 0. The nominator equals zero when s* = 0, 
its derivative equals 
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Expression in the square brackets equals zero when s* = and has derivative 



Clearly, 



dg{s\l3) - 



2/3" e"" (/3- 1) < 0, 



ds* 

and finally g(0,/3) = -1 + 2/?^ - /3^+i < for all /3 G (0,1) because 5(0,1) = and 
•iaiMl = /3^-i[x(l - ;3) + - /3] > 0. Therefore '^^^^ < 0. 

According to Lemma [H Lp can coincide with or only. In case (a), (A) < A 

because lim^^^^ s* = (see l(3T|) ). Function increases as vq decreases (Lemma [3]), but 
$^(wo) = A when vq = d < 13 A. Thus, Vwq S [/9A, A) $^(uo) < A, (iiT + 1) instant jumps 
are never needed and (p = . 

(b) In this case, d > PA according to (a). Since $^(d) = A and $^ is a decreasing 
function (Lemma[3]), $^(uo) > A if -wq € c'] and (^(wo) = *i'^'*'^(^'o), as K jumps are 
not sufficient. Obviously, Lp{d) = <^'^+^(d) = (3 A. Now 

^{pA) = <i>^+i(/3A) = $^+i(d) _ < ci>^+i(d) + ^d-pA)^d 

J PA dvo 

according to Lemma O and statement (a) is proved. 

In case (/?), ^^{vo) < A, hence (^(wo) = ^^{vq). We know that $^((i) = A. Using 
Lemma m we conclude that 

$^(A) = $^(d) + ^^^^dvo >A-iA-d)^d. 



Corollary 4 Theorem[I\ and Corollary hold for undipped cycles. 

Proof. (See FiglHl) Under conditions (a) of Lemma lU (p has a stable stationary point 
V2 coincident with that of (Note that <^^{A) < A, so that V2 € [PA,A).) 
Consider case (b) of Lemma [H 

If vo & [(3A, d] then tp = $^+^ is a contraction defined on this interval; so that the 
statement follows. 

If Vq e {d, A), <p has a stable stationary point V2 coincident with that of . (Note that 
<i>^{A) < A, hence d<A = $^(d), so that V2 G (d, A).) 
Corollary [1] is obvious. H 
Critical cycles. 

Remind that a cycle is called critical if min^ y{s) = From (|25l28l29p it is clear that the 
minimum is attained at 

^o = i^rz^' (39) 

where si solves l(29|) . 
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Lemma 5 Suppose, an undipped k-cycle exists, 
(a) y(so) increases with A. 

(h) For cycles of order k = \, 3e > 35 > 0; ^^^jj^ > e as soon as A > + S. 
Consequently j/(so) ^ oo as A —> oo. 



Proof, (a) After rewriting (|29l) in the form 

1-/3' 
we obtain: 

dsi _ _ (l-e-"i)2(l-/3'=) 



" (40) 
The denominator has derivative (wrt si > 0) 



and hence increases starting from when si — 0. Therefore ^ > 0. 
Since 

y(so) = il + A- vo)e-''° + sq - I + vq - q = vq + sa - q (41) 

we conclude that 



dy{so) 


f dvo 


l_ dso\ 




dA 


\dsi 


dsij 





l-/?*^ si(l-e-''i) J dA 



(b) Note that the denominator in (|40l) is a bounded function of si. Thus 3e > 35i > 0: 
^ > e as soon as si > Si, or, equivalently, as soon as A > + S, where 6 > exists 

because si monotonically increases with A. Remember that Hm . _j_ si = 0. H 

Lemma 6 Suppose, all parameters, apart from b, are fixed. 

(a) A critical cycle of order k exists (for some positive value of b) if and only if 

Y^<<l<<ll (42) 

where ql is given by US]) . The corresponding value of b equals 6o,fe, see ^Bj). 

(b) The boundary ql satisfies inequalities 

^ fc^i <<ll< Al. (43) 



1-/3' 

(In case k = 1, ql = +00.) 
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Proof, (a) Necessity. Let k > 1 and suppose a critical cycle of order k exists. Then, if we 
increase b up to 6* = AJ! — q, this k-cyde (equipped with an asterisk) must remain undipped 
(Lemma [5|): 

y*{s*) = v* + s*-q>0 (44) 
(see (jUJ ), ie q < vX + Sq. Here Sq = In — (see [39]) ), s* solves ((29|) under At and hence 

1 — e 1 

coincides with defined by lfT6|) : is defied by l(28|l . Therefore, + ■^o — 9fc- 
Obviously, system of equations I|28l29l39p and 

Vo + Sq - q = 

(see l|4ip ) must be compatible, ie equation 

'.(«)^%^ + ln^-, = (45) 



must have a positive solution. One can easily check that h increases to infinity with si 

j4gr - q- Hence q > 



starting from lims^^o h{si) — j^r^k ~ Q- Hence q > jzrpi 

In case k — 1 we put ql = +oo, so that (021) transforms to g > and the proof of the 



latter inequality remains unchanged. 

Before proving sufficiency, we firstly prove part (b). 
(b) Let A: > 1; 

^ = '^^ - = T^(-'= + 1) + 1- T^-t;: - T^Tji^ 

= In — a{Tk + 1) - a(Tfc + 1)7 + Tfc7, 

1 — e ^fc 

where a = '"^ i-^fc , 7 = ■ Using l|16|) . the last expression can be rewritten as 

h=\ + In +11- I 7 + Tkl- 

For k > 1 one can easily check that 7 > j:^; therefore, since 1 — e""^*" — r^e"'^'' > 0, 

, ^ , Tk . Tfe 1 - e"^'' - Tfee"'^'' a 

h> 1 h m ■ 



1 - e-^" 1 - e-^" 1 - e-^" 1 - 2a 

(46) 

~ 1 - e-'^'' 1 - e-^fc 1 - e-^fc 3 - Se-'^fc - Tfe - r^e-'^fc ' 

(We have used lfT6|) to express a in terms of r^.) 
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During the proof of Lemma [2ljc), we established that Tk increases with a G (0,1/2), 
starting from when a — 0. Hence Tk G (0, r), where r is the single positive solution to 
equation 

(l-e-^)(l + i(r + l)) = r. 

(The solvability was established in the Proof of Lemma [2jc).) 

Now the righthand side of (|46| is non-negative if G (0, t) . This statement was accuratly 
checked numerically; the analytical proof is problematic. 

The second inequality, to be verified, is obvious: 



^l-^k = + 1) + In - = ^^[^ + "(^^ + ' + < 



(a) Sufficiency. Suppose inequalities (|42l) hold. Then for b G [0, Al — q] (undipped) 
/c-cycles exist according to Theorem [H see FigfTSl (Remember that Al ^ ql — +oo.) Note 
that, in case k > 1, q < Al due to (b). In this case, for b = b* = Al — q, 

y*{4) ^Vo+s*„-q^ql-q>0 

(see l(44|) ) and this particular cycle is really undipped. In case k — 1, according to Lemma 
mb), ?/(so) > for sufficiently large b. Now, if b decreases then the minimal value of y over 
a cycle decreases (Lemmata) ) and, being continuous, becomes zero, since y(so) < for 
the undipped A:-cycle corrresponding to 6 = 0. 

To calculate the critical value of b, note that equation (|l5)) has a single positive solution 
si. Now, if we take 

, si , /3^(si + l) . si .Sl 
b = — + 1 - 9 = — - - In ■ 



1 - e^"! 1-/3'' 1 - e^^i 1 - e-"'! 

then, according to I|28l29p . the corrresponding cycle will be critical. (One can easily see that 
b > 0.) It remains to notice that equation is identical with ifT^ 



Corollary 5 Let N be defined by Then critical cycles of orders k < N cannot exist. 

Proof. According to ifTTI) . q < j^:^, ii n < N. The statement follows from Lemmata). 

■ 

Clipped cycles. 

Proof of Theorem [Tl Let S be the single positive solution to equation 

il + b + S)e-^ = 1. 

Then a continuous trajectory l(25|) starting from (jjq — b, vq = q — S) touches the axis y = 
at a single point, at time moment S. 
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(a) In case S > q — (3 A it is obvious that starting from any point (?/o = b,VQ <E [PA, A)), 
the trajectory never touches the axis y — Q- The statements follow now from Corollary [H 
the mappings ip and (p coincide. 

(b) Suppose that S < q — (3 A and q — S < V , where V {— Vi or V2) is the minimal 
stationary point of the mapping (p (see Lemma|4]and Fig fT4|) . Then, starting from any point 
(2/0 = b^VQ G [(3A,A)), at most Lp(ip{vQ)) is such that the further trajectory never touches 
the axis y — Q- see Lemmas [3] and HI To put it differently, </j"(wo) > g — S* for n > 2. The 
required statements again follow from Corollary 01 The mappings p and (p coincide on the 
domain [q — S,A). 

(c) Suppose that S < q — (3A, q — S > V2, where V2 is the maximal stationary point of the 
mapping p (see Lemma|4]and Fig fT4|) . Then, starting from any point {yo = 6, wo € [Z?^, A)) 
Vn > 2, p^{vo) = (p{(p{vo)) because p{vo), p{p{vo)) < V2 < q — S. Note that in terms of 
Theorem[Tl d < (3{q + b), V2 = p{p{vo)) is different from (smaller than) V2 shown on FigdU 

(d) Suppose that S < q — [3 A, case (b) (Lemma [4]) takes place and Vi < q — S < d (see 
Fig fT4|) . Then, if vq G [f3A,d\, the situation is similar to (c): Vn > 2 (^^"(wo) = p{p{vq))^ 
because p{vq), p{p{vo)) < Vi < q — S < d. 

If Wo S (rf, A) then the trajectory never touches axis y — because \fn p"'{vo) — ip"{vo) > 
d> q — S. The statements follow from Corollary 01 

(e) Suppose that S < q — (3 A, case (b) (Lemma 01) takes place and d < g — 5 < V2 (see 
Fig fT4l) . Then situation is similar to (b). Starting from any point (yo ~ b, vq £ [(3A,A)), at 
most ip{p{vo)) is such that the further trajectory never touches the axis y = 0, the mappings 
p and p coincide on the domain [q — S, A) and the required statements follow from Corollary 



Corollary dl is now obvious. 

Corollary 6 If a clipped k-cycle exists then an undipped k-cycle exists, too. (See ^2^.) 

Proof. As is clear from the proof of TheoremlH < S < q — f3A and p{q—S) = ^''{q — S). 
To put it differently, the domain of is non-empty, so that the corresponding stationary 
point Vi or V2 (Fig fT4l) does exist and defines the undipped /c-cycle. H 

Corollary 7 The order of a clipped cycle can be N or N + I only (see ). 

Proof. Suppose all parameters are fixed, apart from b. For very small values of b, obvi- 
ously, only a clipped A^-cycle is realised. Conditions when a cHpped {N + l)-cycle exists, 
are left till the next subsecion. 

Suppose iV > 1. When we increase 6, fc-cycles with k < N appear: see FigfTSl If b 
is close to Y^pk ~ Q then the fc-cycle has a very short continuous part. From the proof of 
Lemma O we have 



01 ■ 



lim So = and lim y(so) 



b 
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(See (|39|4ip . Therefore, using Lemmata) we conclude that all fc-cycles remain undipped 
indeed. See also Corollary [5l H 

Effects of the router buffer b. 

The goal of this subsection is to justify all the statements of Section HI 
Case A*pf_^^^ < q is trivial: see Fig[l3l Lemmas l5|6| Corollary [7| and its proof. 
Case q < q^^i- According to Lemma [GJ here the {N + l)-cycle appears and becomes 
critical before it extincts at 6 = A*^_^_^ — q. 

Consider the continuous trajectory ((25l) staring from (yo — 0,vo — q): 

( y{r) = e-'' + r -1; 
\ v{r) = q + r. 

Clearly, there is 1 — 1 correspondance between parameters r and b given by equation 

e"'' + r - 1 = 6. (47) 
The {N + l)-cycle cannot be reaHsed if 

P'^iq + r + 1) < y{r) +q^e^'' + r-l + q^b + q. 
Let us study the difference 

A{r)^e-'' + r-l + q-p'^{q + r + l). (48) 

Since ^^7-^ = 1 — — (3^ ^ this difference has a minimum at r = — ln(l— /3^) (corresponding 
to 6 = C, see lfT3|) ) which equals 

q{\ - 13'') - 2/3^ - (1 - /3^) ln(l - /3^) = (1 - /3^)(g - D), 

see (fT2l) . Since the critical {N + l)-cycle exists, we are sure that q < D and the values b 
and b l(20|) are well defined. These equal the minimal and the maximal values providing 
A(r(6)) = 0. Here and below, r(6) is the positive solution to l|47p . Note that the clipped 
(A'' + l)-cycle appears when b = b and becomes critical at 6 = bo^N+i- The value b does not 
play any role because b > &o,Af+i- 

Case q%^i < q < ^jv+i- Here the (A^+l)-cycle cannot be critical (LemmalU). According 
to Lemma m it also cannot be undipped because undipped cycle becomes critical when b 
decreases. Sometimes (A^ + l)-cycles are not reaHsed at all. Firstly, the latter happens if 
D < q. But even if D > q, it can happen that b > A*j^^-^ — q, so that the (A^ + l)-cycle does 
not exist in view of Corollary [6l 

Lemma 7 Suppose < q < ^at+i- 

(a) For a given value of b, the clipped {N + l)-cycle exists iff A(r(6)) < and b < 

(b) A(r(A^^.i - q)) > 0. 

(c) Suppose that D > q. Then b > ^jy+i — q iff C > ^^_f_i — q; b < ^jy+i ^ q iff 
C < ^^+1 - q. 
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Proof, (a) The necessity is obvious: see Corollary [6] and FigfTSl 

Suppose A(r(6)) < and b < A'}^^-^ — q. For the uncHpped {N + l)-cycle, the minimal 
value of y is negative; let us denote the corresponding minimal value of v by v. Then, 
starting from (i/o = 0, uq = <l), the trajectory l(25|) reaches the level y = b, and, after {N + 1) 
instant reductions of v, reaches point {y = b,v < v). After that, the trajectory goes down 
up to the axis y = 0, and the clipped {N + l)-cycle is well defined. 

(b) Value b — A*j^_^_^ — g is the largest buffer size when the undipped (N + l)-cycle 
exists: see Fig[T3l The corresponding minimal value ymin is negative and, starting from 
ivo = ymin, vo = 2/mm + ?) trajcctory (|25l) reaches level y = b a± such value of v that 
j3^v = b + q. Therefore, starting from {yo = 0,vo = q), trajectory l|25p reaches level y = b 
at a smaller value of v, and smaller than {N + 1) reductions of v are needed, meaning that 
A(r(6)) > 0. 

(c) Obviously, b<C<b. Thus the necessity is trivial. The sufficiency follows from (b) 
because A*j^_^^ ^ 9 ^ [k,b]- B 

Corollary 8 In case q%_^^l < q < ^n+h 1 — th^ value of C cannot equal A*j^_^^ — q. 
The proof follows directly from statement (b). Lemma [71 

Corollary 9 Suppose N is fixed. 

(a) For all q G {q'^_^_^, A*j^^^] the value of N remains unchanged, 
(h) IfD> A%^^ then Vg G (g^+i, A^+J C > A%^^ - q. 

(c) Ifqpf+i < D < ^jv+i ^'^s" either C > A*pf_^_^ — D and Vq E [q'^^ijD] C > A'^^-^ — q, 
orC < - D and G (g^^+i, D] C < A*^^^ - q. 

(d) If (3 e (0, 1) varies, equality C = A*j^^^—D can hold only in the area where D < q$f^i. 

Proof, (a) The assertion follows from inequalities 

pN gN-1 
\- - 1n+1 < ^ ^_ pN-1 ' 

see Lemma [2ld) and Lemma [6jb). As usual, ^w-i = +co if = 1. 

(b) Clearly, if q — A*^^-^ = inm{A*j^^^, D} then C > = A*j^_^_^ — q. If we decrease q 
up to the values of C and ^JV+i remain unchanged and situation C — A*j^^^ — q is 
excluded due to Corollary [8l 

(c) The proof is similar to (b): take q = D = minjA^^j^, D} and reduce its value. 

(d) In case C = A'^j^-^ — D and D > q^^i we have a contradiction to (c). H 

One can show that different situations studied in Lemma [7] and Corofiaries [8] and [9] can 
really take place. 

Theorem [2] follows directly from Section [H 

Proof of Proposition [TJ According to definition (fTT)) , S = ^^^^^p-r'^^ , where r solves equa- 
tion 

rjl + P) ^ -r 
1 + 2/3 + /3r ^ ■ 
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The both functions on the left and on the right increase from zero, and t is smaller than 9 
which solves equation -^^2p+)ie = l;ieT<6' = 2/3+1. Now 

rjl+P) ^ . „ -(20+1) ^ ^ (l + 2/3)(l-e-(^/^+i)) 
1 + 2/3 + /3t 1 + /3e-(2/3+i) 

and 

/3 (l + 2/3)(/3 + l)e-(2^+i) 
^>T3^- TT^^^CTO ^ooas/3^1. 

■ 

Proof of Theorem [3l 

First we consider the case b G [0, 6o,i]- this case, the cycle is clipped or critical (see 
Figure [7]). According to Condition (b) of Theorem [2 ii q > the cycle does not have 
multiple jumps for any size of the buffer. Without loss of generality, we assume that the 
zero time moment corresponds to the time moment just after the jump (Point A). Recall 
that we denote the transformed time by s and the original time by t. We denote by Sa the 
transformed time when the system reaches point A, by Sb the transformed time when the 
system reaches point B, and so on. Without loss of generality, we assume that Sa = 0. We 
also use the notation: Sab = Sb ~ Sa = Sb, Sbc = Sc — Sb, and so on. 

From II25D we have 



y{Sc + u)^ycD{u)^e ^ + {u-l), for ue[0,ScD], 

so that 

yiSo) = e"^^" +ScD-l = b. 

We note that v{Sc) — q- Consequently, v{Sd) = 9 + Scd, v{Se) — q + Scd + 1 and 
v{Sa) = P{q + Scd + !)• Again, from l(25|) we have 

y{s) = (1 + g + y{SA) - v{SA))e-' + s-l + v{Sa) - q, 

and 

y{SB) = MS a) + l + q- v{SA)]e~^^^ + [Sab - 1] + v{Sa) - g - 0. 
Thus, we have the following equation for Sab 



B 



[b - (3ScD + (1 - /3)(1 + g)]e-5-^- + Sa 

+/35cD-(l-/3)(l + <z) = 0. 

Now, we can calculate the cycle duration in the original and transformed times. Denote 
these quantities by Tcyde and Scyde, respectively. Note that Scyde — si + 1 (see ((29)) with 
k — 1). From equation v{Se) — v{Sa) + Scyde we obtain 

Scyde = il-l3)iq + ScD + l), 
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and, consequently, 



dt = 







S cycl e 







T 



x{s) 



ds = TS, 



cycle 



m I B 



M Jsa 



y{s)ds 



Sd 



y{s)ds 



Sc 



Next, we calculate the average queue size 



1 



cycle 



mT 



Tcycle Jo 
Sb 



x{t)dt 



T, 



cycle Jo 



x{s) [T + 



x{s) 



y{s)ds 



Sa 



y{s)ds 



Sc 



M \Jsa 



Sb 



y'^{s)ds ■ 



ds 



y\s)ds] +B(T 



Sc 



B 



Now we calculate the average sending rate 



1 



T, 



Using we have 



A = 



cycle Jo 



,{t) 



X{t)dt 



dt 



^ cycle Jo ^cycle Jq 



v{s)ds 



T, 



cycle JO 



{P{q + 1 + Scd) + s)ds 



cycle ^ 



For the calculation of the average goodput we use the following formula: 



9 



T, 



cycle 



JTa \ ^^ / ^ cycle J Sa 



v{s)ds + q + b 



In case b G (&o,i,oo) the cycle is undipped. Consequently, the calculations of the average 
quantities are more straightforward than in the previous case and are based on the knowledge 
of only one parameter S cycle- H 
Proof of Proposition [2j 

If -B — > oo (equivalently, 6 — > oo), then si — > cx) (see equation (1291) 1. According to 
Theorem m we have 

m r f^^ 

q{si + 1) + / y{s)ds + b 



T, 



cycle 



m 



1 + 26 + 2gf- si - (1 + 6 + g 



+ 1) 



1-/3 ' 2 1-/3 

and, consequently, 

^^(2si + 1) - 1 - 26- 2g + si + (1 + 6 + 9 - %±ii)e-«i + ^ 



A = /Lt 



l + 26 + 2g-si-(l + 6 + g 



1-/3 1^ 1-/3 ^ 2(1-/3) *1 
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(2 + ^).si_ 4 1 



l + /3si 



0+, as si oo. 



Proof of Proposition [HI (a) Suppose N is fixed and q = ^ changes, i.e., increases starting 
from TTgjv- Using (fHl) . (fT5|) and omitting for brevity N as the power and the index, we 



obtain: 



1-e- 



(IBq dbo d9 dq ^ 
dm dO dq dm 



{l-e-»)2 



6l(l-e-") 1-/3 



In- 



1. 



We used the impHcit differentiation theorem for Note that ^ — q and express g using 

dBo 



dm 



1-e- 



- hi- 



1 -e- 



In 



- 1 



l-e-o ^ 1-/3 



fl-e- 



l-e-f-ee- 



The second square bracket, f{jS-p), is a monotonous function of 



^ 

(l-e-f) ^ 1-/3 
/3 



(a) If {9 + 1) ^ In iJ^-ii > then /(•) does not decrease and hence 



1 -e 



< — 
dm 1 



-hi- 



- hi- 



-61 



1 -e- 



- 1 - 



1 - e- 

hl-; 



- 1 - /(O) 



1 -e- 



1 



-hi- 



1 - e-" 1 - e- 



(l-e-«) 

7— 1— 7ln7<0, 



A 



because 7 — j— ^ ^ (IjOo) for 9 > and function 7 — 1 — 7ln7 has the maximum which 
is equal to zero at 7 = 1. 
(/3)If 

A-e-^-9e-'^ ,9 
+ r«<0 (49) 



then /(•) decreases and hence 

dBo 



< 



dm 1 



1-e- 



-hi 



1 -e- 



^ - 1 - Urn f{y) 
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e 



1 - e- 

Using (091), we have 



In- 



(1 - e-'^ - ee-<^){e - 1 + e-%{e + 1) 



1 -e 



dBo 

dm 1 — e~ 



' + l)(l-e- 



(1 -e-«)2 



ee- 



6»(l-e-«) 
9e-«)(0-l + e-«)(0 + l) 



3e- 



(l-e-«)2 



(l-e-»)2 

Indeed, consider function g{9) = 3e^^ + d^e^^ 

dg 



< 0, if e > 0. 

2-6^29^ Clearly g(0) = 0; 



^6* 



= 0; 



d[0 + 2e- 



9=0 



de 



= 1 - 2e-^ -29 <0 



because the latter function decreases starting from —1 at 6 — 0. 
Note that 

^-0- as0-O + . 
dm 



(50) 



(b) Obviously, without loss of generality we can put N = 1 and prove that Bq^at increases 
as /3 e (0, 1) decreases. Like previously, we omit A'' as the power and the index. Now again 
using the implicit differentiation theorem we obtain 

dBo _ dbo de 
~dj3 ~ ^~d0 ' dp 



1+9 
'(1-/3)2 



I eil-e-o) 





1-13 A 



< 0. 



Proof of Corollary [31 The first part follows directly from Proposition [31 if we notice that N 



remains unchanged on intervals m g 

"T-i+l- 



and increases by 1 at points 



When m ttin — 0, g approaches and 9n goes to zero (see (15)). According to 



(fTS]). 6o,Ar 0+, hence -Bq^at = mbo^N 0+. Equality (f50|l implies that 



dm 



0-. 
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and ^ 

til i— /a 



In 

Suppose m — > 0+, 5 00, iV = 1, 6*1 ^ cxd. Then — ^ 

according to . Therefore 

Bo.i = mbcAT = — Oi — . 

q \ 9i J 13 

Consider Bo^N{mN-i), ie put q — ^w-i and study equations (fH)) . (fT5|) . Let 9n he the 



positive solution to 



In 



N 



N-1 



N 



then 



Bo,N{mN-i) 



-In- 



7 AT 



1 — e""" 1 — e^''" 



Obviously, liniAr^oo 9n = 0, hence, directly from (|5T|) we obtain 
"In « 



lim 



= - lim — ^ = lim 

2 JV^oo JV^oc 



and finally 



lim mpf-iBQ^pf{mpf^i) = lim 



- In- 



7JV 



1 — e^^" 1 — e^"" 



(51) 



- 1 



N 



1-13 



- 1 



xjirn^l^^) /3Vr)\lim (1-/3^-1)2 = ^ 



2(1-/3) 



/32(A*T)2 = i(l-;3)2(^T)2. 
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